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ABSTRACT 



FUSE and HST/STIS spectra of the dwarf nova WZ Sge, obtained during and 
following the early superoutburst of July 2001 over a time span of 20 months, 
monitor changes in the components of the system during its different phases. The 
synthetic spectral fits to the data indicate a cooling in response to the outburst 
of about 12,000K, from ^ 28, OOO/sT down to ^ 16, OOOiT. The cooling time scale 
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T (of the white dwarf temperature excess) is of the order of 100 days in the 
early phase of the coohng period, and increases to r ~ 850 days toward the end 
of the second year following the outburst. 

In the present work, we numerically model the accretional heating and sub- 
sequent cooling of the accreting white dwarf in WZ Sge. The best compressional 
heating model fit is obtained for a I.2M0 white dwarf accreting at a rate of 
9 X lO~^M0yr~^ for 52 days. However, if one assumes a lower mass accretion 
rate or a lower white dwarf mass, then compressional heating alone cannot ac- 
count for the observed temperature decline, and other sources of heating have to 
be included to increase the temperature of the model to the observed value. We 
quantitatively check the effect of boundary layer irradiation as such an additional 
source. 

Subject headings: Cataclysmic variables - stars: individual (WZ Sge) - white 
dwarfs. 



1. Introduction 

WZ Sge is a well-studied dwarf nova (DN, a sub-class of Cataclysmic Variables - CVs) 
with extreme properties. It has the largest outburst amplitude, shortest orbital period, 
longest outburst recurrence time, lowest mass Roche-lobe filling secondary, and lowest ac- 
cretion rate of any class of DNe (Howell et al. 1999, 2002). In addition it is the brightest 
DN and arguably the closest CV, with a distance of only ^ 43 pc (Thorstensen 2001, private 
communication). WZ Sge was observed to go into outburst in 1913, in 1946 and in 1978 
and it was consequently assumed to have an outburst period of about 33 years. The 23 July 
2001 outburst, first reported by T. Ohshima (see Ishioka et al. (2001)), was therefore 10 
years earlier than expected. 

Following the nomenclature of Patterson et al. (2002), the optical fight curve (Figure 
1) of the system during the July 2001 outburst exhibits a "plateau" phase (this phase is 
really a slow decline phase), which lasted for about 24 days. During this period the sys- 
tem brightness underwent a steady decline with a rate of about 0.1 mag/day, falling from 
My f« 8.2 to My ^ 10.7 in about 24 days, a sign that the accretion was taking place at a 
slowly decreasing rate. It was then followed by a sharp drop (itself lasting a few days) of 
the visual magnitude from m„ 11 down to rriy ~ 13, where it stayed for about 3 days: 
the "dip" . During the dip, accretion had either stopped completely or dropped considerably. 
Between day 29 and day 52 of the outburst, the system then underwent 12 successive "echo" 
outbursts: the " rebrightening" phase. On the 53rd day, the system started to cool without 
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any other noticeable outburst event: the "coohng" phase (in Figure 1 the hght curve is 
shown for t < 75 days only, since for t > 53 days My is a monotonously decreasing function 
of time, as the star cools). 

The 2001 outburst hght curve of WZ Sge is remarkably similar to its 1978 outburst hght 

curve with practically the same initial decline rate of O.lmag/day (Bailey 1979). However, 
in 1978 the plateau and rebrightening phases lasted about 30 days each against 25 days 
each in the 2001 outburst (i.e. the outburst in 1978 lasted about 10 days longer than in 
2001). This may be related to the fact that in 2001 WZ Sge erupted after only 23 years 
of quiescence, while in 1978 it erupted after 33 years of quiescence. The 2001 outburst was 
not as strong as in 1978. Prom the lUE archive we found that the continuum level of the 
spectrum has a flux density Fx (estimated around A = 1, 700 A ) of about Fx 2 x 10"^^ 
ergs s~^ cm~^ during the early phase of the 1978 outburst, or about 3 times larger than 
in the 2001 outburst at the same epoch and wavelength F\ !v 7 x 10~^^ ergs s~^ cm^^ A^^. 
During the cooling phase, the flux density F\ was also larger in 1978 than in 2001 at the 
same epoch (Slevinsky et al. 1999). The 1946 eruption, however, was much different, with 
no apparent rebrightening phase, while the plateau phase had small amplitude variations 
during the whole duration of the outburst. 

The main purpose of the present work is to try to account for the observed sequence 
of temperatures using a numerical model for the accretional heating and subsequent cooling 
of the accreting white dwarf. In the next section we address the issue of the accuracy of 
the temperature determination during the outburst and cooling periods, together with an 
overview of the estimates of the WD temperature and accretion rate. We present the code 
we use to carry out the simulations of the heating and coohng of the WD in section 3. The 
results are presented in section 4. In the last section we discuss the possible origin of the 
observed high temperatures and slow cooling of the white dwarf. 

2. The Accretion Rate and White Dwarf Temperature Determination 

In the present work we consider the observations of WZ Sge carried out in the Far 
Ultraviolet for which the temperature of the accreting white dwarf and/or the mass accretion 
rate of the system were assessed. These observations monitor changes in the FUV component 
of the system during its different phases, over a time span of 20 months, and reveal a cooling 
of the WD in response to the outburst of about 12,000K. Namely, we consider 4 FUSE and 
11 HST/STIS observations as follows: one FUSE observation obtained during the plateau 
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phase, one HST observation obtained during the dip, two HST and one FUSE observations 
obtained during the rebrightening phase and 2 FUSE with 8 HST observations obtained 
during the coohng phase. A total of 15 observations together with their references are 
recapitulated in Table 1. 

In order to estimate the temperature one generates a grid of theoretical spectra for 
different values of the surface temperature, gravity and composition, using the synthetic 
spectra generator codes TLUSTY and SYNSPEC (Hubeny 1988; Hubcny ct al. 1994; 
Hubeny & Lanz 1995). One then masks regions of the observed spectrum that are not 
characteristic of a WD atmosphere (such arc emission lines for example), and a fitting 
technique is used to find the best fit model. The mass accretion rate is estimated in a similar 
manner using different options in the synthetic spectra generator codes TLUSTY (the older 
version is known as TLDISK) and SYSNPEC. 

Mass Accretion Rates: the following mass accretion rates were estimated in the 
references mentioned in Table 1 and should be considered only as an order of estimate. 
From Table 1, we see that the mass accretion rate estimate during the plateau phase is 
rather modest, of the order of 1 — 3 x 10~^ MQ/year (Long et al. 2003a) on day 7 into 
outburst (Table 1). During the dip, the white dwarf might have a temperature of about 
25, OOOK (Knigge et al. 2002) on day 27. This is more than 10, OOOK above its quiescence 
temperature which is ^ 14, 500K. During the rebrightening phase, the mass accretion rate 
seems to be somewhat smaller than during the plateau and only peaks at 1 x lO^^MQyr^-^ 
(Long et al. 2003a) on day 46. We put these values into the table, though we feel that 
they should be considered only as rough estimates rather than actual values and they refiect 
our inability to assess M (and T^^) during this epoch. This is because during the different 
phases of the outburst, additional components (such as the accretion disk and other obscuring 
material ejected during outburst) contaminate and veil (Long et al. 2003b) the spectrum of 
the white dwarf. However, during the cooling phase, one expects to 'see' mainly the white 
dwarf with little or no contribution from the disk and/or other (possibly masking) additional 
components. During this epoch the white dwarf is exposed and its temperature decreases. 

Because of that, we decide to use for our modeling the values of the white dwarf tem- 
perature obtained during the cooling phase only. However, discrepancies of up to 5, OOOK 
exists between the temperature estimates of Sion et al. (2003a) and Long et al. (2003b), for 
the observation carried out on day 50, when the white dwarf was revealed, but apparently 
also partially masked or/and possibly with an accretion disk component. The rest of the 
observations in the cooling phase are consisten within 2, OOOK. 

Temperatures: we elaborate here a little more on the temperature, using the three 
different approaches (denoted a, b and c in Table 2). 
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These were all observations using HST/STIS. In (a) we mask the N V region of the spec- 
trum when needed (together with less affected regions varying from spectrum to spectrum ) 
and use the latest version of the synthetic spectra generator codes TLUSTY and SYNSPEC 
(Hubeny 1988; Hubeny et al. 1994; Hubeny & Lanz 1995); in (b) the spectral fitting 
technique uses the same masking regions for all the spectra together with the new version 
of the code; and (c) is taken from Long et al. (2003b), which uses a third and different 
masking technique together with an earher version of the code. In Table 2, we also hst the 
flux integrated over the wavelength over the entire spectral range of STIS, to the power 1/4, 
for all the epochs, which is proportional to the effective temperature (though each integrated 
flux can be relatively over- or under- estimated). 

In Figure 2 we draw the temperatures hsted in Table 2 as follows: is represented 
by stars; is represented by squares; Tc is represented by plus signs; and the temperature 
estimated through the flux is represented by circles, where it has been arbitrarily scaled so 
as to fit the last data point for Tf,. The two triangles represent the FUSE data points. The 
September 2001 datum has temperature estimates ranging from ^ 27, OOOK up to ~ 32, OOOK. 
For that particular epoch the spectrum was taken between the last two echo outbursts and 
obviously the white dwarf was not the only component and in addition might have been 
partially masked. For that reason we decide not to take into consideration the data points 
that were not known accurately, or which were taken during periods when the white dwarf 
might still have been veiled. In our numerical modeling we will try to model the cooling 
temperature as listed in Table 2, with some reservations about the September 2001 entry. 

The cooling of the white dwarf is initially pretty rapid and then slows down significantly. 
If we assume an exponential cooling law T = Tinf + AT x e"*/^ (with Ti^f = 14,500K), 
equivalent to AT{t) — AT x e"*/"^, where AT{t) is the excess temperature of the WD at 
time t, then we find that r is not constant, i.e. the cooling is not exponential. In Tables 1 and 
2, t is counted from the begining of the cooling period, day 52 into outburst, however, it is 
usually agreed to count the cooling time from the beginning of the cooling phase. Therefore, 
we set here t ^ t — 52 days, and wc denote the cooling time at time t = n days by r". We 
find ~ (70 +io) days, ^ (190 i30) days, r^^^ ^ (325 l40) days, and r^^^ « (850 t280) 
days. 

It is interesting to note that the observed cooling following the 2001 July outburst is 
faster than the cooling observed after the 1978 outburst. Slevinsky et al. (1999) found a WD 
temperature of 20,500K, 220 days after the outburst against ~ 19, OOOK at the same epoch 
for the 2001 outburst (e.g. observing epoch of April 2002, day 266). And ~ 500 days after 
the 1978 outburst the temperature was still around 17,500K while the present observations 
seem to indicate a lower temperature (fa 17, OOOK) 466 days after the 2001 outburst. 
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3. Numerical Modeling 



In order to model the accretional heating and subsequent coohng of the WD in WZ Sge, 
we use a ID quasi-static evolutionary code without hydrodynamics (quasi-static assumption). 
It is an updated version of the quasi-static stellar evolution code of Sion (1995). The 
code includes time variable accretion, OPAL opacities and boundary layer irradiation which 
indirectly accounts for the stellar rotational velocity (see Eq.l below). Stellar rotation, 
however, is not included anywhere else in the code. In this code, the white dwarf is computed 
all the way down to the core, in the region well below the nuclear biirning region. Initial 
models are constructed by the fitting point method and the resulting initial model down to a 
depth where p = lO^g cm'^ is stored as input for the evolution code. For a given white dwarf 
mass and an initial effective surface temperature (Tgjj), we chose R^d to be the theoretical 
equilibrium radius. The models built in this manner are initially in equilibrium, and if no 
matter is accreted onto the WD surface, the WD parameters do not undergo any change on 
the time scale studied here. All other details of the code can be found in Sion (1995) and 
references therein. 

Numerical simulations are carried out by switching on accretion for the duration of 
a dwarf nova outburst interval and then shutting it off to follow the cooling of the white 
dwarf. In this way the effects of compressional heating and boundary layer irradiation can 
be assessed quantitatively. The matter is assumed to accrete 'softly' with the same entropy 
as the white dwarf outer layers. From theoretical considerations, one expects the accreting 
matter, as it transits through the boundary layer region, to increase its temperature (since 
the boundary layer temperature is high: Tbl » T^). In addition in the boundary layer 
region, the radial velocity is large and the accreting matter can also heat up by shock. For 
these reasons the 'soft' accretion assumption is actually not justified, since this advected 
energy is non-negligible. However, in our treatment of the boundary layer irradiation, we 
take into account part of this energy by assuming that a fraction of the energy liberated in 
the boundary layer (Lbl) is absorbed by the outer layer of the star: it is included as a source 
in the energy equation in the outer layer of the star. 

The treatment of the boundary layer irradiation is done as follows. The energy liberated 
in the boundary layer is given by (Kluzniak 1987): 




(1) 



where 



acc 



GM^M 
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is the total accretion energy, G is the gravitational constant, M* is the mass of the star, 
is the radius of the star, M is the mass accretion rate, is the angular rotation rate of the 
star and VtxiR*) is the Keplerian angular velocity at one stellar radius. Equation (1) can be 
used as long as the disk is geometrically thin and optically thick and extends to the stellar 
surface. In the present case one expects the stellar rotation rate to be the rotation rate of 
the WD: fi* = fi^^. 

We assumed that only a fraction of the boundary layer luminosity is irradiating the star, 
namely: 




A value a — 1 means that half of the BL luminosity is lost into space while the other half is 

absorbed by the star. Assuming a value a — 0.5 means that only 25% of the BL luminosity 
is absorbed by the star. Here, we choose a = 0.5, which is the value used in the work of 
Shaviv & Starrfield (1987), and about half the estimated value of Regev & Shara (1989) 
who used a ~ 0.2Lacc- We then assume different rotation rates 

For a non rotating star 77 = and the BL luminosity is exactly half the accretion energy 
Lbl — Lacc/'^, while for a star rotating near break-up r) — 1 and Lbl — 0. In this work we 
took values ranging from rj = 1 (for no boundary layer irradiation, in order to check only the 
effect of the compressional heating) down to rj = 0.05 (when the star is slowly rotating and 
boundary layer irradiation is the main source of heating). Clearly a smaller value of a will 
require a smaller value of rj in order to keep the same amount of BL irradiation in a specific 
model, and vice versa. However, the compressional heating results obtained in this work are 
not at all affected by the value of a used in the simulations. 

In the simulations it is assumed that the accretion and heating of the white dwarf is 
uniform rather than being restricted to the equatorial region. In addition, the transfer of 
angular momentum (by shear mixing) into the white dwarf is neglected. 

4. Results 

As previously remarked in section 2, during the outburst, as accretion takes place at 
a high rate, the star's photospheric emission is overwhelmed by the emission of the hot 
components (mainly the inner disk), which makes it difficult to assess the exact temperature 
of the star and its rotation rate ftwd- However, on day 53 the system is found in a low 
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state and starts to cool. During that time the accretion rate has probably dropped to its 
quiescence level. 

Numerical simulations (Godon & Sion 2002) have shown that the temperature increase, 
due to BL irradiation, is sustained only during accretion, and when the accretion is turned 
off, the star rapidly radiates away the BL energy absorbed in its outermost layer. However, 

the temperature increase due to compressional heating takes place deeper in the outer layers 
of the star and it takes many days (months) for the star to cool. Therefore, we assume that in 
the cooling phase the observed elevated temperature of the star is due to the compressional 
heating it has endured during the outburst phase alone. 

Since BL irradiation and compressional heating take place at different depths in the outer 
layers of the star, and on different time scales, we model their effect on the temperature of 
the star in separate model runs. 

First we model the compressional heating of the star alone ocuring during the outburst 
phase, therefore setting r) = Qyjd/^K = 1 in the simulations. The exact mass accretion rate 
of the system is not known and it is changing with time during the outburst. From the optical 
light curve in Figure 1, it seems very likely that initially the mass accretion rate is very high 
at the onset of the outburst and decreases steadily during the plateau phase. At present, our 
code only simulates a constant accretion rate, but at any time the accretion can be turned off 
to model the cooling. Consequently, we model the outburst phase (plateau + rebrightening) 
by accreting at a constant rate for 52 days, after which accretion is turned off and the WD 
cools. From Table 1, we expect the mass accretion rate during the outburst phase to be of 
the order of (say) 3 x lO^^MQyr"^ and the temperature of the star to increase from about 
15, OOOK up to 28, OOOK (or more, see Table 2) around day 53. We run models with different 
white dwarf mass, namely M^^ = 0.7Mq, LOMq, L2M0 and vary the mass accretion rate 
in the range lO^^MQyr"-'^ < M < 10~^ MqJY^^ . We chose the initial WD temperature to be 
around 13, 500 — 15, OOOK. Ideally we are looking for a model that reaches a temperature 
of 28,000ir (or more) after 52 days of steady accretion, with M 3 x lO"~^M0yr~^, and 
with a temperature curve as close as possible to the observations of Tables 1 & 2. In Table 3 
we hst models (1-4) obtained for different values of Myj^ and M. From these models, model 
3 with M^d = 1-2 X Mq and M = 9 x lO-^Mg/year, yields the best fit. Both models 1 
& 2, with a lower white dwarf mass {0.7Mq and 1 Mq respectively), and model 4, with a 
lower mass accretion rate (5 x 10~^MQ/year), yield a lower temperature and a sharper initial 
temperature drop than the observations. In figure 2 we draw the temperature increase of the 
WD due to compressional heating alone, together with the observed values, as a function of 
time for model 3. Models 4 & 2 are drawn in Figures 3 & 4 respectively. 

Clearly, from Figures 3 & 4, we see that compressional heating alone cannot account 
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for the observed temperatures if the mass M^d < 1.2 x Mq and/or the mass accretion rate 
M < lO~^M0/year, because it produces lower temperatures than observed. Since both a 
1.2Mq mass and a lO~®M0/year accretion rate are larger than expected, we decided to 
check the effect of boundary layer irradiation on the temperature of the white dwarf for a 
model with M = I.OMq in addition to the compressional heating. Namely, we decide to 
check how boundary layer irradiation increases the temperature of the model in addition 
to compressional heating. Specifically, we chose model 2 arbitrarily around t—175 day. 
There it has a temperature of only 18,000K (due to compressional heating. Figure 4) while 
the observations indicate that the temperature should be at least 2,000K higher at that 
epoch. We therefore run a 1.0 Mq model (model 5, same mass as model 2) with an initial 
temperature of 18,000K (to simulate day 175 of model 2) and a quiescent mass accretion rate 
(M ?a 10^^^ — lO^^^MQ/year), and we include boundary layer irradiation assuming a rotation 
rate f^^^^ = 0.2^2^^ (corresponding to the observed rotational velocity of 1200km/s). Namely, 
we wish to check if the additional heating of the boundary layer irradiation of the accretion 
disk during the cooling phase (quiescent M) can increase the temperature of the model by 
2,000K, and for what value of M (here we model the combine effect of compressional heating 
and boundary layer irradiation separately, model 2 and model 5 respectively). 

We find, from model 5 in Table 3, that one can increase the temperature of model 2 at 
t=175 days by 2,000K by means of boundary layer irradiation if the mass accretion rate is 
at least M Ki2x IQ-^^Mg/year. 

5. Discussion and Conclusion 

In the present numerical exploration we first assumed that all the accretional heating 
of the white dwarf was due to compressional heating, and we found that the mass of the 
accreting white dwarf needed in the simulations of the WZ Sge sup erout burst to fit the 
observations was rather large: 1.2Mq. The average mass accretion rate of the outburst 
model, 9 x lO^^Mgyr"^ , assessed from the compressional heating simulations, was itself larger 
by a factor of f« 5, than the average value determined from the spectral fits to the observations 
which was ~ 1.7 x lO^'^Afoyr^^ [25 days of accretion at a rate M = 3 x lO~^M0yr~^, followed 
by 25 days of intermittent accretion at an average rate of ~ 5 x lO~^°M0yr~^]. 

When we assumed a smaller WD mass of only one solar mass, we found that even with 
a high mass accretion rate of lO~^M0/year the model could not account for the observed 
temperature. For that model, we simulated a second source of heating (in addition to 
compressional heating): boundary layer irradiation at a quiescence mass accretion rate, 
corresponding to the mass accretion rate during the cooling period. The mass accretion rate 
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we obtained (2 x lO^^^Mo/yr; Table 3) to fit the observations is in agreement with a recent 
analysis of an X-ray observation of WZ Sge. Hasenkopf & Eracleous (2002) presented results 
of a uniform analysis of all the ASCA X-ray observations of non-magnetic CVs. For WZ Sge 
they estimated an X-ray luminosity of about Lx ~ 2.7 x 10^° erg s~^ in the range 0.5-10 keV, 
assuming a distance of 69 pc. Rescaling this value to the better estimate of 43 pc (and also 
to be consistent with the distance assumed in this work) leads to Lx ~ 7 x 10^° erg s~^ for 
the X-ray Luminosity. Hasenkopf & Eracleous (2002) found that optically thin boundary 
layer models (Popham 1999) provide the best description of the data. Since the disk is 
radiating in both +z and -z directions, the X-ray luminosity is at least half the boundary 
layer luminosity, namely: Lx — Lbl/'^i and if the star does partially mask the inner disk 
and boundary layer (since i = 78°), then one has Lx < Lbl/'2- For this reason we assume 
Lbl ~ 3 X Lx ~ 2 X 10^^ erg s^^, and using Kluzniak (1987)'s relation for Lbl, sji angular 
velocity of 1200km s^^ and a 1.0 solar mass, we find the quiescence mass accretion rate from 
the X-ray observation to be M = 6 x lO~^^M0/yr, about 3 times smaller than our estimates. 

Another possibility, that we are unable to assess quantitatively, is the slow release of 
rotational kinetic energy from the outer layer of the star (Sparks ct al. 1993). If the outer 
layer of the star (equatorial belt) has been spun up during the outburst, then rotational 
kinetic energy can be release during the cooling phase as this fast rotating layer spins down. 
This effect could account for an additional source of heating of the white dwarf ocuring 
during the cooling phase. 

In this work we have shown that compressional heating alone can account for the cooling 
curve of WZ Sge following the July 2001 superoutburst only if the WD in WZ Sge is massive 
{Myjd = 1.2Mq) and the outburst accretion rate is large (^ lO'^Mgyr"^). If this is not the 
case, then compressional heating alone is not enough to account for the observed decrease in 
the WD temperature, and we suggest boundary layer irradiation from a quiescent accretion 
disk and the slow release of rotational kinetic energy from a fast rotating accretion belt 
as possible additional sources of heating of the white dwarf to account for the observed 
temperatures. 

We are thankful to Knox Long for providing us with temperature estimates from his 
work before publication. We are grateful to the VSNET collaboration for providing us 
with the data shown on Figure 1. BTG acknowledges support from a PPARC Advanced 
Fellowship. P. Godon is particularly thankful to C. Hasenkopf for providing us with a copy 
of her paper on the X-ray analysis and to Thierry Lanz for kindly providing us with the 
latest version of the Tlusty and Synspec codes. S. Starrfield acknowledges partial support 
from NSF and NASA Grants to ASU. 
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Table 1. FUSE & STIS observations of WZ Sge 



obs. 


date 


day 


instrument 




M 


phase 


Reference^ 


number 








< 1, ooox > 


< M^/yr > 






1 


30 Jul 01 


7 


FUSE 


- 


l-3e-09 


plateau 


[1] 


2 


19 Aug 01 


27 


STIS 


25.0 


- 


dip 


[2] 


3 


22 Aug 01 


30 


STIS 


- 


5e-10 


rebrightening^ 


[2] 


4 


07 Sep 01 


46 


FUSE 


42.0 




rebrightening^ 


[1] 


4 


07 Sep 01 


46 


FUSE 




le-09 


rebrigntenmg 


ri 1 
[1] 


5 


11 Sep 01 


50 


STIS 


Tabic 2 


3e-10 


rebrightening^ 


[3] 


6 


29 Sep 01 


68 


FUSE 


25.0 




cooling 


[1] 


7 


10 Oct 01 


79 


STIS 


Table 2 




cooling 


[3] 


8 


07 Nov 01 


107 


FUSE 


23.0 




cooling 


[1] 


9 


10 Nov 01 


110 


STIS 


Table 2 




cooling 


[3] 


10 


11 Dec 01 


141 


STIS 


Table 2 




cooling 


[3] 


11 


15 Apr 02 


266 


STIS 


Table 2 




cooling 


[3] 


12 


05 Jun 02 


317 


STIS 


Table 2 




cooling 


[3] 


13 


27 Aug 02 


400 


STIS 


Table 2 




cooling 


[3] 


14 


01 Nov 02 


466 


STIS 


Table 2 




cooling 


[3] 


15 


23 Mar 03 


608 


STIS 


Table 2 




cooling 


[3] 



^References: [1] Long et al. 2003a; [2] Knigge et al. 2002; [3] see Table 2. 

^The observations on day 30 and day 46 coincide with the first and 10th echo outburst of the 
rebrightening phase respectively. The observation made on day 50 was made during a relatively 'low' 
state between the 11th and the 12th outburst in the rebrightening phase. The observation on day 
46 was modeled twice: once with a white dwarf atmosphere only and once with an accretion disk 
only. The observation on day 50 was modeled once with a white dwarf atmosphere only and once as 
a combination of an accretion disk and a white dwarf atmosphere. 
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Table 2. STIS temperature estimates of WZ Sge 



obs. 


date 


day 


^1 

a 


-'-b 


c 


(Flux)V4 


number 






< looo/r > 


< looo/r > 


< looo/r > 




5 


11 Sep 01 


50 


31.9 


27.0 


28.2 


3.8x10-^ 


7 


10 Oct 01 


79 


25.2 


23.6 


23.4 


3.3x10-=^ 


9 


10 Nov 01 


110 


23.7 


22.4 


22.1 


3.1x10-3 


10 


11 Dec 01 


141 


22.6 


21.6 


20.7 


2.9x10-3 


11 


15 Apr 02 


266 


19.5 


18.8 


18.1 


2.5x10-3 


12 


05 Jun 02 


317 


18.8 


18.0 


17.4 


2.5x10-3 


13 


27 Aug 02 


400 


17.8 


17.4 


16.7 


2.4x10-3 


14 


01 Nov 02 


466 


17.5 


17.0 


16.3 


2.3x10-3 


15 


23 Mar 03 


608 


17.2 


16.6 


15.9 


2.2x10-3 



^Ta and Ti, were estimated in this work (see text) while Tc is from Long et al. 
2003b. 
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Table 3. Accretional heating models of WZ Sge 



Model 
number 


< Mq> 


LogRyjd 
< cm > 


< 1, ooox > 


M 

< Mq/vt > 


rpmax{2) 
-'-wd 

< 1, OOOX > 


^wd 
<^K > 


lengtli(^) 
< days > 


1 


0.7 


8.95 


15.0 


le-08 


21.0 


1.0 


52 


2 


1.0 


8.78 


15.0 


le-08 


25.0 


1.0 


52 


3 


1.2 


8.60 


14.5 


9e-09 


34.0 


1.0 


52 


4 


1.2 


8.60 


14.5 


5e-09 


27.5 


1.0 


52 


5 


1.0 


8.78 


18.0 


2e-ll 


20.0 


0.2 





^The initial temperature of the white dwarf. 

^The maximum temperature reached by the model during the run. 

^Models for which the effect of boundary layer irradiation is taken into account {^wd < 1) are 
run for a few days only. 
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Figures Caption 

Figure 1: The VSNET optical light curve of the 2001 outburst of WZ Sge. The first 6 
observations are marked by arrows. The four different phases of WZ Sge have been marked 
for clarity. 

Figure 2: Modeling the heating and cooling of WZ Sge. The temperature (in Kelvin) of 
the white dwarf is drawn as a function of time (in days) since the start of the outburst (July 
23, 2001) The sohd hne represents the compressional heating modeling, for 1.2 solar mass 
white dwarf with an initial temperature of 14,500K, accreting at a rate of 9 x lO~^M0/year 
for 52 days. The stars, squares and plus signs denote the temperatures listed in Table 2 Ta, 
Tjj and respectively. The triangles are the FUSE data points and the circles represent the 
temperature estimates using the flux values (see text). 

Figure 3: Same as Figure 2, but here the solid line represent a model with a mass 
accretion rate of 5 x 10~^Mo/year and an initial temperature of 15,000K (model 4 in Table 
3). 

Figure 4: Same as Figures 2 and 3. The solid line represents model 2 in Table 3. 



-17- 



- 




1 ' 1 ' 1 


1 


- 






CO — ^ . «J. 






- 








- 






"'it',; 


ng 








.-;«'** 


o 










o 




- 






o 


- 


- 

- 

- 




- ■ • 




- 

- 

- 


- 

- 




^ ' "I 




- 
- 








CTi 










f — 








.> ' 
































_c 










cn 










' ^ 










_Q 








"' 






- 




.. 

"w 






— 




n — > „; 




- 








CL 


















# 






- 








- 












- 








- 






















- 






a 


- 


















o 










Q_ 
















#■ 




















1 


, 1 , 


1 





o 

CD 



O 
LO 



o 

^ A 
cn 

o 



CD 



O 

CN 



O 



o 



00 O CN CD 

AIAI 



-18- 




